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Abstract 

Wilson quark action, with tadpole improved clover term added, is studied on coarse 
anisotropic lattices. The bare velocity of light parameter in this action is determined 
non-perturbatively using the pseudo-scalar and vector meson dispersion relations for 
various values of the gauge coupling f3 and bare quark mass parameter k. 
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1 Introduction 



X 



It has become clear that anisotropic, coarse lattices and improved lattice ac- 
tions are the ideal candidates for lattice QCD calculations on small computers 
[1,2,3,4]. They are particularly advantageous for heavy objects like the glueballs, 
one meson states with nonzero three momenta and multi-meson states with or 
without three momenta. The gauge action employed is the tadpole improved 
gluonic action on asymmetric lattices: 
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where Pij is the usual plaquette variable and Rij is the 2x1 Wilson loop on 
the lattice. The parameter Ug, which we take to be the 4-th root of the average 
spatial plaquette value, incorporates the so-called tadpole improvement and ^ 
designates the (bare) aspect ratio of the anisotropic lattice. 

Using this action, glueball and hadron spectrum has been studied within the 
quenched approximation [3,4,5,6,7,8,9]. The next step would be to extend the 
calculation to more complex physical quantities. However, in order to utilize the 
improved fermionic action [10,11,12] on anisotropic lattices, some parameters in 
the action have to be determined, either perturbatively or non-perturbatively, in 
order to obtain as much improvement as possible. In this letter, we will discuss 
the tuning of these parameters in a quenched calculation using clover- improved 
Wilson fermions on anisotropic lattices. This paper is organized in the following 
manner. In Section 2, the clover-improved Wilson fermion action on anisotropic 
lattices is introduced. In Section 3, we focus on the non-perturbative tuning of 
the bare velocity of light, using the dispersion relations of pseudo-scalar and 
vector mesons. The optimal values of the bare velocity of light are obtained for 
various values of f3 and k on different lattices. In Section 4, we will conclude 
with some general remarks. 



2 Improved Wilson Fermions on Asymmetric Lattices 



Consider a finite four- dimensional lattice with lattice spacing a^ along the /i 
direction for /i = 0, 1, 2, 3. For definiteness, we also denote ag = at and aj = 
tts for i = 1,2, 3. We will use ^ = ttg/at to denote the bare aspect ratio of 
the asymmetric lattice. The un-improved Dirac- Wilson operator on asymmetric 
lattices is given by: 

Do = Iy: ^m[(Vm + Kh, - a,v,v;] , (2) 



where V^ and V* are the usual lattice forward and backward derivatives. The 
parameter z/^ enters the Dirac- Wilson operator since the hyper-cubic symmetry is 
broken. In particular, we have z/q = 1 and Ui = u which we call the "bare" velocity 
of light. The clover-improved Dirac- Wilson operator is obtained by adding a 
Scheikholeslami-Wohlert, or clover, term to the usual Wilson fermion action [13]. 
To implement the tadpole improvement, one replaces each spatial link Ui{x) by 
Ui{x)/us while keeping the temporal link unchanged. In quenched calculations, 
one usually needs to calculate the quark propagators at various valance quark 
masses. This amounts to different values of itlq or k for the same gauge field 
configuration. In this case, it is convenient to use the following fermion matrix: 



A =S 



''^^xv ^xv^ ~r ^xy 

3 

l/(2/tmax) + Pt X! ^Oi^Oi + Ps(o'l2^12 + 0-23-^23 + 0-31-^31; 
i=l 

- H ^M [(1 - 7M)f^M(^)^=«+M,J/ + (1 + 7M)f^i(a; - P)4-M J ' (3) 

where the coefficients are given by: 

r]i = u/{2us) , r]o = i/2 , a = 1/{2k) - l/{lKmax) , 

Pt = csw{l+0/i^u^s) , Ps = csw/i2ut) . (4) 

In this notation, the bare quark mass dependence is singled out into the param- 
eter a and the matrix A remains unchanged if the bare quark mass is varied. 
Therefore, one could utilize the shifted structure of the matrix Ai to solve for 
quark propagators at various values of itlq (or equivalently k,) at the cost of solv- 
ing only one value of /t = Kmax, using the so-called Multi-Mass Minimal Residual 
(M^R) algorithm [14,15,16]. 

3 Tuning of the parameters 



Tuning the bare parameters of the fermion matrix to obtain the best improve- 
ment is a rather delicate issue [11]. In a quenched calculation, the matter is 
simplified since the anisotropy ^ can be determined from the pure gauge sector. 
In fact, with the tadpole improvement, the renormalization effects of the aspect 
ratio can be ignored [3,4]. The parameter csw is the hardest to be tuned. In 
principle, one could utilize the Schrodinger representation of the QCD path in- 
tegral to determine this parameter non-perturbatively. This has been done for 
the clover- improved Wilson fermion action on symmetric lattices [17,18]. In this 
work, the tadpole improved tree-level value is used. In what follows, we will 
focus on the tuning of the parameter z/. Since this parameter has rather weak 
dependence on the quark mass on anisotropic lattices [11], it can be determined 
non-perturbatively by the pseudo-scalar meson dispersion relation for a rather 
heavy pion. 

Results of this tuning process for various values of (3, together with the corre- 
sponding simulation parameters are given in Table 1 to Table 5. The value of 
f3 in our simulations ranges between 1.7 and 3.0, which roughly corresponds to 
physical spatial lattice spacing a^ in the range 0.12 ~ 0.4/m. The aspect ratios 
of these lattices are fixed at ^o = 5. The spatial sizes of the lattices range from 
L = 4 to L = 8, which in physical units lie in the range 0.5 ~ 3.2fm, depending 
on the value of (3. For simplicity, only the results for 6^40 lattices that are close 
to the optimal value of z/ are shown. 



Table 1 

Simulation parameters for /? = 1.7 lattices. Input aspect ratio parameter ^ is fixed 
to be 5 throughout this paper. The approximate spatial lattice spacing a^ in physical 
units as obtained from [4] is also shown. Results for pion mass, rho mass in lattice 
units, the slope of the pseudo-scalar and vector meson dispersion relation Z-,^ and Zp 
(see Eq. (5) for the definition of them), and the corresponding value of v and k are 
also listed. 



/3 = 1.7, as ~ 0.39/m, 6^40 lattices 


V = 0.90, 288 configurations 


V = 0.92, 192 configurations 


K, 


"T-TrOi mpttt Zt, Zp 


K 


m.„at rupat Zt, Zp 


.0585 
.0575 
.0565 
.0555 
.0545 


0.3029(8) 0.416(2) 0.97(2) 1.01(5) 
0.3803(7) 0.477(1) 1.06(1) 1.08(4) 
0.4472(5) 0.537(1) 1.09(2) 1.15(4) 
0.5159(4) 0.597(1) 1.18(1) 1.06(4) 
0.5827(4) 0.657(1) 1.23(1) 1.08(4) 


.059 
.058 
.057 
.056 
.055 


0.230(1) 0.360(5) 0.97(2) 1.04(6) 
0.315(1) 0.427(2) 1.02(3) 1.12(5) 
0.3916(9) 0.491(4) 1.11(2) 1.19(5) 
0.4625(8) 0.552(3) 1.17(2) 1.10(7) 
0.5319(7) 0.612(2) 1.22(2) 1.15(5) 



Table 2 

Same as Table 1 except for /? 



2.2. 



/3 = 2.2, as ~ 0.27 fm, 6^40 lattices 


u = 0.90, 160 configurations 


z^ = 0.95, 208 configurations 


K 


m^af nipat Z^, Zp 


K 


m^fli mpttt Zt, Zp 


.061 
.060 
.059 
.058 
.057 


0.198(1) 0.291(3) 0.89(4) 0.93(6) 
0.268(1) 0.341(2) 0.96(3) 1.00(5) 
0.3361(7) 0.397(2) 1.00(3) 1.07(4) 
0.3994(5) 0.452(2) 1.07(2) 1.16(4) 
0.4616(5) 0.509(1) 1.19(2) 1.18(5) 


.060 
.059 
.058 
.057 
.055 


0.203(1) 0.294(2) 0.99(4) 0.98(6) 
0.276(1) 0.351(2) 1.04(3) 1.07(5) 
0.3443(6) 0.406(2) 1.10(2) 1.17(4) 
0.4091(5) 0.463(1) 1.18(2) 1.25(4) 
0.4728(4) 0.520(1) 1.26(2) 1.34(3) 



Table 3 

Same as Table 1 except for /? 



2.4. 



/3 = 2.4, as ~ 0.22/m, 6'MO lattices 


V = 0.90, 192 configurations 


V = 0.92, 192 configurations 


K 


m^Ci nipat Z^, Zp 


K 


m^at nipat 2"^ Zp 


.061 
.060 
.059 
.058 
.057 


0.207(1) 0.265(3) 0.96(3) 1.05(5) 
0.2713(9) 0.320(2) 1.08(2) 1.13(4) 
0.3351(8) 0.377(1) 1.16(2) 1.19(3) 
0.3971(7) 0.434(1) 1.22(2) 1.26(3) 
0.4579(6) 0.491(1) 1.28(2) 1.27(3) 


.0605 
.0595 
.0585 
.0575 
.0565 


0.204(1) 0.265(3) 0.99(3) 1.10(6) 
0.273(1) 0.319(2) 1.08(3) 1.17(4) 
0.338(1) 0.376(2) 1.16(3) 1.22(4) 
0.3973(8) 0.433(1) 1.26(3) 1.28(4) 
0.4587(6) 0.491(1) 1.31(2) 1.31(4) 



Table 4 

Same as Table 1 except for (3 



2.6. 



f3 = 2.6, as ~ 0.19/m, 6^40 lattices 


V = 0.85, 192 configurations 


V = 0.88, 192 configurations 


K 


niT^at rUpat Z^, Zp 


K 


m^at nipat 2"^ Zp 


.062 
.061 
.060 
.059 
.058 


0.206(1) 0.243(2) 0.99(3) 1.08(4) 
0.2687(8) 0.301(1) 1.07(2) 1.16(3) 
0.3288(6) 0.358(1) 1.16(2) 1.22(3) 
0.3886(8) 0.415(1) 1.20(2) 1.26(3) 
0.4468(7) 0.471(1) 1.25(2) 1.30(2) 


.0615 
.0605 
.0595 
.0585 
.0575 


0.199(1) 0.236(2) 1.02(3) 1.09(5) 
0.2628(8) 0.295(2) 1.11(2) 1.21(3) 
0.3243(6) 0.354(1) 1.18(2) 1.22(3) 
0.3852(8) 0.411(1) 1.25(2) 1.32(3) 
0.4446(7) 0.469(1) 1.31(2) 1.36(3) 



Quark propagators with various three momenta are obtained using the M^R 
algorithm with wall sources. In this work, only the meson energy for zero and 
the following three non-zero lattice momenta are measured, namely (100), (110) 



Table 5 

Same as Table 1 except for /? 



3.0. 



/3 = 3.0, as ~ 0.12/m, 6^40 lattices 


u = 0.85, 192 configurations 


V = 0.88, 192 configurations 


K 


m^Oi mpttt Zt, Zp 


K 


m^ot mpOt Zt, Zp 


.0635 
.0625 
.0615 
.0605 
.0595 


0.149(2) 0.152(2) 1.00(3) 0.90(4) 
0.199(1) 0.206(2) 1.02(2) 1.09(3) 
0.253(1) 0.262(1) 1.10(2) 1.16(2) 
0.3086(9) 0.318(1) 1.17(2) 1.22(2) 
0.3644(8) 0.373(1) 1.16(2) 1.27(2) 


.062 
.061 
.060 
.059 
.058 


0.191(1) 0.202(2) 1.05(2) 1.09(3) 
0.246(1) 0.258(1) 1.14(2) 1.16(3) 
0.302(1) 0.314(1) 1.21(2) 1.24(2) 
0.3593(9) 0.371(1) 1.20(2) 1.29(2) 
0.4161(8) 0.427(1) 1.26(2) 1.34(2) 



and (111) with Dirichlet boundary conditions applied in the temporal direction. 
Then the pseudo-scalar correlation functions and the vector meson correlation 
functions are constructed, from which the single pion energy and the single rho 
energy are extracted in lattice units. Due to the advantage of the anisotropic 
lattices, the meson energy levels can be measured with good accuracy. Errors 
are obtained using the standard jack-knife method. The mass of the pion and 
rho, are also listed in Table 3. The single meson energy levels are fitted according 
to the expected dispersion relation: 



ul{\.)=ml + ZX , colik) 



m^ + Zpk' 



(5) 



In Fig. 1, the dispersion relations of the pseudo-scalar (the lower plot) and the 
vector meson (the upper plot) are shown for j3 = 2.2, u = 0.95 and five different 
values of k on 6^40 lattices. The straight lines represent the fits using the first 
three lowest momentum points (including the zero-momentum point) according 
to Eq. (5). It is seen that the above dispersion relations fit the simulation data 
well for lattice momenta that are not too large. Fitting qualities are similar 
for other parameter sets. The results of the slope obtained from the fitting are 
also included in Table 3. By tuning the bare velocity of light parameter u, we 
obtain a value of Z^ which is consistent with unity. The same tuning procedure 
is also applied to smaller (4^40) and larger (8^40) lattices. The outcome of the 
procedure is quite similar. No significant dependence of the parameter u on the 
volume of the lattice has been observed. 



4 Conclusions 



In this letter, we outline the method on the non-perturbative determination of 
the bare velocity of light parameter u in the anisotropic Wilson quark action Q. 
The optimal values of u are obtained for various values of (3 and k. These values 
can then be utilized in a quenched calculation using the improved anisotropic 



^ We have just noticed that, in a recent work, similar tuning process has been studied 
by H. Matsufuru et al. [19]. 



^ = 2.2 v=0.95 6^40 Lattice 




Fig. 1. The dispersion relation for the pseudo-scalar (the lower plot) and vector meson 
(the upper plot) obtained on 6^40 lattices with f3 = 2.2 at five different pion masses. 
The bare velocity of light is taken to be 0.95. It is seen that the measured dispersion 
relations agree with the free dispersion relation at small lattice momenta. 

Wilson action. Some of these calculations are in progress and we hope to report 
the results soon. 
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